Flat Chern Band in a Two-Dimensional Organometallic Framework 
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By combining exotic band dispersion with nontrivial band topology, an interesting type of band 
structure, namely the flat chern band (FCB), has recently been proposed to spawn high-temperature 
fractional quantum hall states. Despite the proposal of several theoretical lattice models, however, 
it remains a doubt whether such a "romance of flatland" could exist in a real material. Here, 
we present a first-principles design of a two-dimensional (2D) Indium-Phenylene Organometallic 
Framework (IPOF) that realizes a nearly FCB right around the Fermi level by combining lattice 
geometry, spin-orbit coupling and ferromagnetism. An effective four-band model is constructed to 
reproduce the first-principles results. Our design in addition provides a general strategy to synthesize 
topologically nontrivial materials in virtue of organic chemistry and nanotechnology. 

PACS numbers: 73.43.Cd, 73.61. Ph 



Whenever an unconventional band structure is brought 
out from a conceptual model into a real-world material, 
a wide range of theoretical advances and technological 
innovations will be triggered. The well-known examples 
are graphene [J and topological insulators , which 
are featured with exotic dispersive bands and nontrivial 
topological bands, respectively. By combining these two 
features, another interesting type of band, namely the 
flat Chern band (FCB) has recently been proposed 0-0] , 
which is dispersionless and characterized by a nonzero 
Chern number. A well-dehned FCB requires its band 
width smaller than both the energy gap between the FCB 
and other bands and the interaction energy scale 0-0] ■ 
Since the kinetic energy is strongly quenched, carriers in 
the FCB experience strong Coulomb interaction in addi- 
tion to the topological frustration that in together spawn 
unprecedented topological strongly-correlated electronic 
states pj-[loj. Fundamentally different from the narrow 
bands commonly existing in heavy fermion compounds, 
which has a trivial Chern number, the outcome of a FCB 
requires a delicate balance of lattice hopping, SOC and 
ferromagnetism. Because of the stringent criteria, no 
real material to date has been experimentally observed 
to contain a FCB. 

Advances on synthetic chemistry and nanotechnol- 
ogy have shown the potential in producing complex 
lattices [ill, E3- Recent experiments using substrate- 
mediated self-assembly have successfully fabricated 2D 
organometallic frameworks with different lattice symme- 
try. For example, Shi et al. [l3[ have synthesized two 2D 
frameworks by co-depositing tripyridyl molecules with Fc 
and Cu atoms on Au (111) surface, which form Kagome 
and triangular lattices, respectively. These covalent or- 
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ganic frameworks are found to exhibit remarkable ther- 
mal stability. In this Letter, we present a first-principles 
design to realize the FCB in a 2D Indium-Phenylene 
Organometallic Framework (IPOF). Density functional 
theory (DFT) calculations [14| show that in this unique 
IPOF structure, a nearly flat band appears around the 
Fermi level characterized by a nontrivial Z2 topological 
number. Upon p-type doping, the spin degeneracy of the 
flat band is spontaneously lifted; the ferromagnetic flat 
band is then found to have a nontrivial Chern number. 

Figure la shows the atomic structure of IPOF, whose 
key feature is to bind p-oribtal heavy elements (In) with 
organic ligands (para-phenylenes) into a hexagonal lat- 
tice. Each hexagonal unit cell contains three phenylenes 
and two In atoms. As a common feature of group-Ill el- 
ements, the In atoms naturally bond to three phenylenes 
in a planar triangular geometry. Such a bond configura- 
tion is identical to that in triphenyl- indium In^CeH^^, 
a common indium compound fla . IT(| . A similar or- 
ganic framework assembled from 1,4-benzeneddiboronic 
acid on Ag (111) surface has been fabricated in exper- 
iment, with only In atoms replaced by boroxine rings 
[PI Ha |. The stability of IPOF has been examined 
by first-principles lattice relaxation and phonon calcu- 
lations; details are discussed in Sec. II of the Supple- 
mentary Material (SM). The proposed synthesis process 
of IPOF is described in Sec. Ill of SM. 

To examine the electronic structure of IPOF, we first 
purposely exclude SOC from the calculation. The result- 
ing electronic band structure, wavefunction and atomic- 
orbital projected density of states (DOS) are summarized 
in Fig. lc. The band structure presents a nonmagnetic 
insulator picture. Electronic bands at the band edge ex- 
clusively come from the p-orbitals of C and In atoms. 
Specifically, at the valence edge, there are ten bands 
within 2eV from the valence band maximum, which can 
be divided into two groups. Around E—Ep = leV, there 
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FIG. 1: (Color online) (a) The atomic structure of IPOF. 
(b) the first Brillouin zone and special k-points. (c) (Left) 
Band structure without SOC from DFT (blue solid curves) 
and model Hamiltonian (eq. (8) with t\ = 0.7eV, A = 0, 
M = (red dashed curves). Insets are the wavefunction iso- 
surfaces of two states denoted by arrows. (Right) Atomic 
orbital projected density of states. 



are six p z bands, which can be attributed to the highest 
7r-electron molecular levels from the three benzene rings, 
as shown by both the wavefunction and the DOS. The 
remaining four bands exhibit the in-plane p xy features, 
which arise from hopping among the cr-bonds between In 
and C atoms (red shaded ellipses in Fig. la). The top 
and bottom bands in the four p xy bands are nearly flat 
in the whole Brillouin zone with a narrow bandwidth of 
lOmeV (without SOC). The middle two dispersive bands 
form a Dirac cone at the K points similar to graphene. 
The flat bands touch the dispersive bands at the T point. 
In contrast to the benzene molecular levels, the wavefunc- 
tion of the p xy flat bands surprisingly has a distribution 
across the whole lattice, i.e. they are not localized states 
(insets of Fig. lc). This is the first indication that these 
flat bands are unique and nontrivial. 

Next, we include SOC in the calculation, and the re- 
sults are shown in Fig. 2a. Comparing Fig. 2a with Fig. 
lc, the most significant difference is that several degen- 
erate points of the p xy bands, e.g. the T and K points, 
are split. Consequently, the two flat bands become sep- 
arated from the dispersive bands. The direct (A^ 2 r ) and 
indirect gap (A* 2 d ) between the top flat band and the 
nearest dispersive band (Fig. 2c) are A^ 2 r = 90meV and 



^■Ind = SOmeV. The separation between the top and 
bottom flat bands ( A 14 in Fig. 2a) is lAeV. Another 
effect of the SOC splitting is that the band width {W) of 
the flat bands increases from 10 to 60meV. Note that all 
the bands are still spin degenerate because of the time- 
reversal and inversion symmetry. The topology of the 
flat band is characterized by a nontrivial Z2 toplogical 
number as discussed in Sec. IV of SM. 

Partially-filled flat band is unstable due to the large 
DOS at the Fermi level. It has been rigorously proved 
that even an arbitrarily small Coulomb interaction will 
drive the system into a ferromagnetic ground state at 
specific filling factors [l9|, [2(| • To obtain some insights 
within the DFT formalism, we dope the system by manu- 
ally reducing the number of valence electrons in the unit 
cell, while maintaining the charge neutrality with a com- 
pensating homogeneous background charge. This makes 
the top flat band partially filled. Calculations reveal a 
spontaneous spin polarization perpendicular to the 2D 
plane. Figs. 2b and 2c shows the band structure of this 
ferromagnetic ground state under doping. The spin-up 
and spin-down bands are separated apart, with the Fermi 
level shifting below the topmost spin-polarized flat band. 
This ferromagnetic ground state is 5meV lower than the 
spin-unpolarized state in total energy, which can be fur- 
ther stabilized by applying an external Zeeman field. The 
spin splitting, U ~ lOOmeV (Fig. 2b), represents the 
strength of the on-site Coulomb interaction. We thus 
estimate the electron-electron interaction of the p xy elec- 
trons in IPOF to be on the same order. For comparison, 
the spin splitting of the p z electrons (benzene molecular 
levels) is one order of magnitude smaller. 

We will now examine the topology of the topmost spin- 
polarized flat band by directly calculating its Chern num- 
ber based on its DFT wavefunctions. The Chern number 
is defined as: 



c=^-[ d 2 kF 12 (k), 
2tt« Jbz 



(1) 



where the Berry curvature i 7 i2(k) is given by: 

f 12 = dMk) - d 2 A 1 (k) 

A. L = (u k \di\u k ) 



(2) 



|uk) is the Bloch function of the flat band and the 
derivative di stands for djd^. The distribution of Berry 
curvature is shown in Fig. 2e, exhibiting an interesting 
ring pattern. The integration of the Berry curvature in 
the whole Brillouin zone gives c = 1. Therefore, we con- 
clude that this top fiat band is a FCB. 

A direct manifestation of the Chern number is the 
number of chiral edge modes circulating around the 
boundary. To check the edge property, we have calcu- 
lated the momentum-resolved edge DOS of a semi-infinite 
IPOF 14]. As shown in Fig. 2d, one single edge band 
emerges in the gap (A 12 ) and unidirectionally connects 
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FIG. 2: (Color online) (a) Band structure with SOC. (b) Band 
structure with SOC when doping one hole into the unit cell. 
The solid (blue) curves are DFT results. The dashed (red) 
curves are from model Hamiltonian eq. (8) with parame- 
ters (a) ti = 0.7eV, X = 0.05eV, M = 0; (b) ti = 0.7eV, 
X = 0.05eV , M — O.leV (c) Zoomed-in band plot around the 
Fermi level, (d) Momentum-resolved edge density of states 
of a semi-infinite IPOF with SOC and doping. The bright- 
ness is proportional to the magnitude of the density of states, 
(e) Distribution of the Berry curvature Fi2(k) of the top flat 
band in the reciprocal space. Left is the DFT result; right is 
the model result using the same parameters as in (b). 



the top fiat band with the dispersive band below it, signi- 
fying the typical feature of a chiral mode and indicating 
the Chern number c = 1. 

We have further constructed an effective Hamiltonian 
to understand the origin of the FCB in IPOF. Consid- 
ering that the fiat band comes from hopping among the 
a bonds, we use these cr-bond orbitals as basis to con- 
struct the effective Hamiltonian. There are three a bonds 
around each In atoms as marked by the red elliptical 
shades in Fig. la. We first symmetrize the basis orbitals 
by considering the nearest-neighbor (NN) hopping (to) 



I I 
H = t Q \ 10 1 

1 1 



(3) 



The eigenstates of H consist of a doublet (ed) and a 
singlet (e s ): 
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The symmetrized eigenstates transform like vectors: 
^ di and v E'd2 within the x-y plane; \I/ S perpendicular to 
the x-y plane. As long as \E' S is far away from the dou- 
blet, it can be reasonably neglected. The problem is then 
reduced to a four-dimensional subspace spanned by two 
sets of ^di and 1 $>d2 associated with the two In atoms in 
the unit cell. Considering the next-NN (NNN) hopping 
(ti) via the benzene rings, we can write out a four-band 
hopping Hamiltonian for each spin: 
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V yy = 4(1 + e ikai + 4e' ik ' a2 ) ; a 1)2 is the lattice vector 
as shown in Fig. la. Without loss of generality, we have 
set the on-site energy of "F^i and ^ d 2 as zero. The eigen- 
values of Hhop give four energy bands (the red dashed 
curves in Fig. lc), which agree well with the DFT re- 
sults. Hhop effectively describes a honeycomb lattice with 
two in-plane vector-like orbitals at each site, i.e. the p xy - 
orbital counterpart of graphene. As discussed in Refs. 
pll |22| . this kind of hopping Hamiltonian guarantees 
two completely flat bands, because their Wannier func- 
tions forbid any electron "leakage" to outside through a 
destructive interference. 

With regard to SOC, we start from the original atomic 
form: Hsoc = XL • a. By realizing that ^di,d2,s trans- 
forms like three vectors, we can write out its second quan- 
tization form: 

H soc = X(ic dH Cd2i ~ «4it c <* 2 t + c !t Cd H - credit 

+icl t c d2 i - ic^c^t + h.c), (6) 

in which cj CT creates an electron with spin a z on the 'Fj 
orbital. To reduce H SO c into the four band (^1^2) 
subspace, we keep only the leading terms: 



(7) 
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which do not include coupling between different spin 
components, so the spin- up and spin-down spaces get au- 
tomatically separated. For each spin subspace, the SOC 
can also be written as a 4 x 4 matrix: 
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where a z = ±1 is the spin eigenvalue. This SOC repre- 
sents an imaginary on-site coupling, independent of mo- 
mentum. It pins an additional phase to the electrons 
when they hop between the two on-site orbitals and re- 
sult in nontrivial topology when electrons travel in the 
lattice via Hy lov . 

The spontaneous magnetization under doping is a 
many-body effect [lj|[2(|. Within the single-electron pic- 
ture, we reproduce the ferromagnetism by simply adding 
a Zeeman term. Finally, the effective Hamiltonian to de- 
scribe the single-electron property of IPOF can be writ- 
ten as: 



ff e //,<7„=±i(*) = H hop (k) + H 
( V xx V xy \ 
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There are three parameters in the model: the NNN 
hopping parameter t±, the spin-orbit coupling strength A 
, and the spontaneous magnetization M. Note that the 
NN hopping term to does not explicitly enter the four- 
band model that renders the flat band more accessible 
in IPOF than in other model lattices 0-01 • The higher- 
order hopping terms, e.g. the next-NNN hopping will af- 
fect the flatness, but they are already negligible in IPOF, 
smaller than 10 meV according to the DFT results. By 
fitting the three parameters to the DFT results, the dis- 
persion of the four p xy bands can be well reproduced, as 
shown in Fig. lc, 2a and 2b. To ensure that this simple 
model also reproduces the wavefunction property, we in 



addition calculate the Berry curvature using the model 
eigenstates; the result is also in good agreement with the 
DFT result (Fig. 2e). The middle two dispersive bands 
from the model are found to have c = 0; the bottom flat 
band has c = 1 ; thus the total Chern number of the four 
band subspace is zero. These justifications confirm the 
validity of the effective Hamiltonian, which will be useful 
for future investigations going beyond the DFT formal- 
ism. 

In summary, we reiterate in Tab. I several key energy 
scales associated with the FCB in IPOF. Further studies 
including many-body effects are necessary to determine 
the ground state and low-energy excitation of the car- 
riers in the FCB. The emergence of fractional quantum 
hall state requires further enlarging the energy gap A 12 
and reducing the bandwidth W to satisfy the condition 
A > U > W [1, 0, H- To achieve this intriguing state, 
the organic nature of IPOF provides the flexibility to tune 
the parameters in kinds of chemical ways, e.g. by func- 
tionalizing the benzenes with different chemical groups, 
replacing the benzenes with other organic ligands, or us- 
ing different metal atoms. Furthermore, similar design 
principles can also be applied to realize other exotic band 
structures. For example, by simply replacing In with Bi 
and Mn, one makes new organic topological insulator and 
magnetic topological insulator, respectively (23[. 



TABLE I: Energy scales associated with the FCB in IPOF 



Property 


Symbol 


Value 


Ref. 


Band width 


W 


QQmeV 


Fig. 2c 


Spin splitting 


V 


lOOmeV 


Fig.2b 


Debye temp. 




300meV 


Fig.Sl 


Energy gap 


A 12 


90meV 


Fig.2c 




A 12 


30meV 


Fig.2c 




A 14 


lAeV 


Fig.2a 
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